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ABSTRACT 

An XMM-Newton observation of the bright QSO PG1211+143 in 2001 revealed a 
blue-shifted absorption line spectrum indicative of a high velocity radial outflow of 
highly ionised gas. Unless highly collimated, the outflow mass rate was shown to be 
comparable to the accretion rate, with mechanical energy a significant fraction of 
the bolomctric luminosity. Analysis of the full XMM-Newton data set now allows the 
wider effects of that energetic outflow to be explored. We find that absorption and re- 
emission of the primary continuum flux in the ionised outflow, together with a second, 
less strongly absorbed, continuum component can explain the strong 'soft excess' in 
PG1211+143 without the extreme velocity 'smearing' in conflict with observed absorp- 
tion line widths. Previously unpublished data from a second XMM-Newton observation 
of PG1211-I-143 is shown to be consistent with the new spectral model, flnding that 
the additional continuum component dominates the spectral variability. We speculate 
that this variable continuum component is powered by the high velocity outflow. 

Key v^rords: galaxies: active - galaxies: Seyfert: quasars: general - galaxies: individ- 
ual: PG1211-I-143 - X-ray: galaxies 



1 INTRODUCTION 

, The luminous QSO PG1211-I-143 has been observed twice by 
XMM-Newton, in 2001 and 2004. Analysis of the EPIC and 
RGS spectra from the first observation provided evidence 
for a highly ionised outfiow with a velocity of ~0.09±0.01c 
(Pounds et al. 2003; hereafter P03), though this high veloc- 
ity was contested following a careful ion-by-ion modelling 
of the RGS data (Kaspi and Behar 2006). However, a re- 
examination of the 2001 EPIC data (Pounds and Page 2006) 
resolved additional absorption lines in the EPIC MOS spec- 
trum, strengthening the initial claim, and yielding a pre- 
ferred velocity of v~0.14c. Confirmation of the high velocity 
outflow is important since, as noted in P03, the mechanical 
energy in the flow, if not highly collimated, is a signiflcant 
fraction of the bolometric luminosity of PG1211-I-143 and 
could be typical of AGN accreting near the Eddington rate 
(King and Pounds 2003), while also providing an example 
of the feedback required by the linked growth of SMBH in 
AGN wfth their host galaxy (King 2004). 

P03 also found a strong 'soft excess' over a 2-10 keV 
power law fit of photon index r~1.8, a typical value for 
type 1 AGN, and a broad Fe K emission line. Subsequently, 
the 2001 EPIC spectrum of PG1211-I-143 has been used by 
Gierlinski and Done (2004), and more recently by Schurch 



and Done (2006), to demonstrate how strong absorption of 
the intrinsic X-ray continuum in a 'velocity-smeared', high 
column, ionised gas could provide an alternative explanation 
(to Comptonisation) for the strong soft excess generally seen 
in type 1 AGN (Wilkes and Elvis 1987, Turner and Pounds 
1989). A more general study by Chevallier et al. (2006) also 
considered an ionised reflection origin of the soft excess in 
AGN (e.g. Crummy et al. 2006) and concluded that absorp- 
tion was the more likely cause of a strong soft excess (as in 
PG1211-H143). 

In the present paper we use the full XMM-Newton data 
set from both observations of PG1211-I-143 to re-examine 
the broad-band X-ray spectrum of PG1211-I-143, to better 
understand the structure and dynamics of the outflow, and 
and thereby to assess the conjecture that the massive and 
energetic outfiow substantially defines the emerging spec- 
trum. 

We assume a redshift for PG1211-M43 of z = 0.0809 
(Marziani et al. 1996). 

2 OBSERVATION AND DATA REDUCTION 

PG121H-143 was observed by XMM-Newton on 2001 June 
15 for -^53 ks, and again on 2004 June 21 for ~57 ks. In 
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this paper we mainly use the high signal-to-noise data from 
the EPIC cameras (Striider et al. 2001; Turner et al. 2001), 
but with a crucial constraint on high resolution spectral fea- 
tures provided by the simultaneous RGS data (den Herder 
et al. 2001). All X-ray data were first screened with the 
XMM SAS v6.5 software and events selected correspond- 
ing to patterns 0-4 (single and double pixel events) for the 
pn camera and patterns 0-12 for the MOS cameras. We ex- 
tracted EPIC source counts within a circular region of 45" 
radius defined around the centroid position of PG1211-I-143, 
with the background being taken from a similar region, off- 
set from but close to the target source. After removal of 
data during periods of high background the effective expo- 
sures for the pn camera were ~49.5 ks in 2001 and ~25.2 ks 
in 2004, with corresponding exposures of ~103 ks and ~81 
ks for the combined MOS cameras. Individual spectra were 
then binned to a minimum of 20 counts per bin, to facilitate 
use of the minimalisation technique in spectral fitting. 

Spectral fitting was based on the Xspec package (Ar- 
naud 1996), version 11.3. All spectral fits include absorp- 
tion due to the line-of-sight Galactic column of Nh ~ 
2.85 X 10^°cm"^ (Murphy et al. 1996) and errors are quoted 
at the 90% confidence level (Ax^ = 2.7 for one interesting 
parameter) . 



3 VISUAL EXAMINATION OF THE EPIC 
DATA 

Figure 1 displays the pn spectra of PG1211-(-143 from the 
2001 and 2004 observations in a way that has become con- 
ventional, by first fitting a power law over the 2-10 keV band 
and then extrapolating this fit to lower energies, revealing 
the 'soft excess' typical of many type 1 AGN. This approach, 
followed by P03 in their initial analysis of the 2001 pn cam- 
era data of PG121 1-1-143, yielded a power law fit over the 
2-10 keV band with a 'canonical' photon index of r~1.79, 
together with an Fe K emission line (apparently with a broad 
red wing) , and a strong soft excess parameterised by a black 
body spectrum of kT~0.1 keV. The same approach for the 
2004 pn data of PG1211-I-143 yields a steeper power law 
(r~1.85) and a less strong 'soft excess' (figure 1). 

An alternative comparison of the 2001 and 2004 spec- 
tra is shown in figure 2, which plots the unmodelled, 
background-subtracted data. Shown in this way, we see the 
spectral change is actually an increase in flux (from 2001 
to 2004) between ~0.7-2 keV, and a (less obvious) decrease 
at ~0.5-0.7 keV. Interestingly, the 2001 and 2004 spectra 
are essentially identical above ~3 keV. We choose to show 
the MOS data in figure 2 since, while the broad flux profile 
is the same, the higher energy resolution of the MOS data 
shows additional low energy structure in the 2001 spectrum. 

To highlight the broad spectral changes between the 
2001 and 2004 data we then directly compared the two data 
sets, noting, for example, that if absorption is having a con- 
trolling effect on the emerging spectrum, it should show up 
in the way in which individual spectra differ. As variable ab- 
sorption (multiplicative) spectral components are best seen 
in the ratio of data sets, we therefore re-binned the EPIC 
data from 2001 to a minimum of 200 counts per channel, 
in order to show any broad features more clearly, and then 
re-binned the 2004 data to have identical energy channels. 
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Figure 1. Comparison of the X-ray spectra of PG1211+143 dis- 
played in the conventional way, showing the strong soft excess 
sitting above a 'canonical' power law, for the pn data from the 
2001 (black) and 2004 (red) observations. On this view it appears 
that the soft excess in 2004 is weaker but 'hotter' 
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Figure 2. Direct comparison of the background-subtracted MOS 
data from the 2001 (black) and 2004 (red) observations shows the 
spectral change to be due to an increase in flux between ^0.7— 2 
keV, and a less obvious flux decrease at ~0.5-0.7 keV 

Dividing the first data set by the second then gave the spec- 
tral ratio plot of figure 3. The most obvious feature is a 
broad flux deficit over the ~0.7-2 keV band, the shape and 
positioning of which is strongly suggestive of differential ab- 
sorption by ionised gas, with enhanced absorption in the 
2001 spectrum. 

Figure 4 shows the 'difference spectrum', obtained by 
subtracting the 2004 (background-subtracted) data set from 
that of 2001. The data have again been re-binned to a mini- 
mum of 200 counts per bin. While above ~2 keV the differ- 
ence spectrum reflects the very similar hard X-ray fluxes, in 
addition to a lower flux at ~0.7-2 keV we see a significant 
excess peaking at ~0.5 keV in the 2001 data. 

In summary, a direct comparison of the 2001 and 2004 
MOS data (comparison of the pn data shows the same fea- 
tures) suggests a decrease in continuum absorption is the 
main cause of the spectral change in PG1211-I-143 from 2001 
to 2004, together with a decrease in an emission component 
specific to the soft X-ray band. Implicitly, both results sup- 
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Figure 3. The MOS spectral ratio obtained by dividing the 2001 
count rate spectrum by that obtained in 2004. The most obvious 
feature is a broad deficit at ~0.7-2 keV indicative of variable 
photoionised absorption 
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Figure 4. The MOS difference spectrum, obtained by subtract- 
ing the 2004 spectrum from that of 2001, shows the near-identical 
hard X-ray spectra, with a lower flux at ~0.7-2 keV and compa- 
rable higher flux peaking at ~0.5-0.6 keV in the 2001 observation 

port the contention that reprocessing in ionised gas is a pri- 
mary factor in determining the form of the observed broad 
band X-ray spectrum. An important constraint on the na- 
ture of the variable absorption is the similarity in the ion- 
ising flux during the 2001 and 2004 observations indicated 
by the near-identical spectra above ^3 keV, suggesting a 
change in covering factor rather than ionisation parameter. 



4 THE CASE FOR A FRESH APPROACH TO 
MODELLING THE X-RAY SPECTRUM OF 

PG1211-I-143 

Accumulation of high quality, broad-band X-ray spectra 
of AGN, particularly iiomXMM-Newton EPIC, has raised 
doubts over the the physical reality of the soft excess and 
also the ubiquity of the relativistic Fe K emission line sug- 
gested by earlier ASCA observations. In particular, the con- 
ventional interpretation of the soft excess in type 1 AGN, as 
arising by Comptonisation of cooler disc photons, has been 
questioned by Gierlinski and Done (2004). 

Those authors studied 26 AGN from the PG sample ob- 
served with XMM, fitting the X-ray spectra with a double 
Comptonisation model. They found that while Comptonisa- 
tion of the accretion disc emission provided a good statistical 
fit to the soft excess below 2 keV in every case, the temper- 



ature of the cool Comptonising region was remarkably con- 
stant over the sample, despite that sample having a large 
range in mass and luminosity, and hence in disc tempera^ 
ture, and in the ratio of power released in the hot plasma 
to that from the disc. Given these difficulties, Gicrinski and 
Done (2004) proposed that the soft excess was an artefact of 
ionised absorption in a sub-relativistic wind above the inner 
accretion disc. In order to explain the generally smooth form 
of the soft excess the absorbing gas was required to have a 
complex velocity structure. Subsequently, Sclmrcli and Done 
(2006) have explored this 'smeared absorber' model in more 
detail, showing - in particular - that by including re-emission 
from the ionised gas the strong soft excess in PG1211-I-143 
could be reproduced. 

Chovallior ct al. (2006) have also examined alternative 
origins of the soft excess in AGN, including enhanced soft 
X-ray emission from an ionised accretion disk (e.g. Crummy 
et al. 2006). Chevallicr et al. found that ionised absorp- 
tion in a high speed outflow, or an inhomogcncous accretion 
flow could best explain some AGN X-ray spectra, particu- 
larly those with very strong soft excesses (as PG1211-I-143), 
while only weak excesses can be modelled by reflection un- 
less the primary continuum is strongly depressed (Fabian 
et al. 2002). However, they concluded that simple 'smeared' 
absorption models would require 'fine tuning' to constrain 
the depth of the trough near 1 keV and thereby be generally 
applicable. 

A second spectral feature whoso interpretation has boon 
questioned recently is the broad wing seen to the low energy 
side of the Fe K fluorescent emission line in a number of 
bright AGN. While widely considered to be due to Doppler 
and strong gravity broadening of radiation refiected off the 
innermost accretion disc (eg Fabian et al. 2000), an alter- 
native explanation as (again) an artefact of absorption has 
been proposed in specific cases (eg Inoue and Matsumoto 
2003, Kinkhabwala 2003). 

Clarifying the nature of such significant features in the 
X-ray spectra of AGN is important to explore the accre- 
tion process, outfiows of ejected matter, and ultimately the 
properties of the central black hole. Its strong soft excess, 
evidence for a massive and energetic outflow, and high ac- 
cretion ratio make PG1211-I-143 a particularly interesting 
object by which to explore an alternative approach to de- 
convolving the broad band X-ray spectrum. 



5 DEVELOPING A NEW SPECTRAL MODEL 
FOR THE 2001 OBSERVATION OF 
PG1211-I-143 

We begin the fresh approach by noting that the intrin- 
sic continuum over the broad X-ray band now visible to 
XMM-Newton might not be adequately described by a sim- 
ple power law derived by fitting to the data over the ~2-10 
keV band. Recent studies of spectral variability in MCG-6- 
30-15 (Vaughan and Fabian 2004) and 1H0419-577 (Pounds 
et al. 2004), among the most detailed broad-band spectral 
analyses to date, have provided clear evidence of a variable 
power law component of slope significantly steeper than the 
'canonical value' of r~1.9 (Nandra and Pounds 1994). That 
finding is important in the present context since, if absorp- 
tion is to create the impression of a strong 'soft excess', then 
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the underlying continuum must be much steeper than it ap- 
pears above ~2 keV. To allow for that possibility, while re- 
taining a harder 'primary' continuum in the range predicted 
by the disc/corona model (Haardt and Maraschi 1991), we 
include 2 power law components in the new model. For sim- 
plicity we will call these the primary and secondary contin- 
uum components. 

Indication from the flux ratio plot that the main spec- 
tral change between the 2001 and 2004 XMM-Newton ob- 
servations of PG1211-(-143 is due to a change in ionised ab- 
sorption affecting the spectrum above ~0.7 keV, suggests 
the model should include a moderately ionised ('warm') ab- 
sorber, in addition to the highly ionised absorber responsible 
for the narrow, blue-shifted absorption lines (P03, Page and 
Pounds 2006) . Finally, we include components to model the 
re-emission from both ionised absorbers. 

In Xspec terms the prospective model is then 
WA(P01*ABS1*ABS2 + P02*ABS3*ABS4 + EMl/3 + 
EM2/4), where WA represents the line-of-sight Galactic col- 
umn, included in all fits, and POl, P02 are the separate 
power law components. ABSl and 3 represent the high ion- 
isation absorber, differing only in column density on the 2 
continua, while ABS2 and 4 similarly represent the lower 
ionisation absorption. Finally, EMl/3 and EM2/4 model the 
re-emission from the highly ionised and warm absorbers, re- 
spectively. All ABS and EM components are modelled by 
photoionised gas using the XSTAR code (Kallman et al 
1996). We assume a relatively high turbulent velocity of 1000 
km s~^ FWHM in order to allow substantial line opacity be- 
fore saturation. Free parameters in each XSTAR grid are the 
column density and ionisation parameter, with outflow (or 
inflow) velocities included as an adjustment to the appar- 
ent redshift of the absorbing or emitting gas. All abundant 
elements from C to Fe are included with the relative abun- 
dances as a further variable input parameter, but tied for 
all the XSTAR components and limited to within a factor 
2 of solar. Further constraints on the model are to tie the 
ionisation parameters of ABSl and 3, and of ABS2 and 4, 
and also with their respective re-emission spectra. 

We first evaluated the model with the pn and MOS data 
sets, separately, given their differing low energy responses 
and the higher MOS resolution, which potentially is impor- 
tant in resolving features in the soft X-ray band (Pounds 
and Page 2006). 

An excellent statistical fit was obtained for the pn 
data with the component parameters as listed in Table 1. 
The spectral components are shown together with the total 
model and pn data in figure 5. Fitting the model to the 2001 
MOS data yielded similar component parameters, although 
the statistical quality of the fit was less good. The fit pa- 
rameters are listed in line 2 of Table 1 and the data, model 
and separate model components are shown in figure 6. 

Visual examination of figures 5 and 6 clarifies the main 
features of the new model. 

Above ^^3 keV the spectrum is dominated by the 'pri- 
mary' power law, with a photon index r~2.2, at the upper 
end of - but consistent with - the 'accepted' range for type 
1 AGN spectra. At lower energies the secondary power law 
component (r~3.1) is more important, while below ~1 keV 
re-emission, particularly from the 'warm' absorber, becomes 
significant. We note that variability in this component be- 
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Figure 5. The full spectral model (red) fitted to the 2001 pn cam- 
era data of PG1211+143. The model components are the primary 
power law (green), secondary power law (dark blue), re-emission 
from the moderately ionised gas (light blue) and re-emission from 
the highly ionised gas (magenta) 
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Figure 6. The spectral model (red) fitted to the 2001 MOS cam- 
era data of PG1211+143. The model components are primary 
power law (green), additional power law (dark blue), re-emission 
from the moderately ionised gas (light blue) and re-emission from 
the highly ionised gas (magenta) 



tween 2001 and 2004 could explain the low energy peak in 
the difference spectrum of figure 4. 

Strong absorption of the primary power law is seen to be 
responsible for the mid-band (~2-6 keV) spectral curvature, 
with ionised gas column densities of N£f~3 — 5 x lO^'^ cm~^. 
In both pn and MOS spectral fits the model shows a broad 
trough near ~0.8 keV, due to Fe-L absorption. Figure 3 
shows a similar spectral feature in the 2001/2004 ratio plot, 
again indicating lower absorption in the second observation. 

The secondary power law continuum appears to suf- 
fer much less absorption, with a negligible warm absorber 
column and a factor ~2-5 lower column density of highly 
ionised absorber. 

The high velocity, high ionisation absorber produces the 
narrow absorption lines reported in P03 (with revised iden- 
tifications in Pounds and Page 2006), while dilution by the 
steep power law component could explain the weakness of 
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Table 1. Parameters of the model fits to the pn, MOS and RGS data from the 2001 XMM-Newton observation of PG1211+143. Power 
law indices Fi and r2 refer to the primary and secondary continuum components, respectively. Both high and moderate ionisation 
absorbers aifect each continuum component, with equivalent hydrogen column densities in units of 10^^ cm~^ and ionisation parameters 
in erg cm s~^. Effective redshifts and implied velocities from the ionised gas absorption and re-emission spectra are given in the text 
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narrow absorption lines in the RGS data (P03, Kaspi and 
Behar 2006). 

5.1 Velocity structure 

The velocity structure of the photoionised gas components 
is of particularly interest and relies on identifying signifi- 
cant spectral structure in the data. Our approach in P03 
and in Pounds and Page (2006) was based on the visual 
identification of narrow absorption lines, backed up by mod- 
elling with XSTAR. While that works well where the aim, 
as there, was to determine the velocity of a strong outfiow, 
here we have to rely principally on modelling to interpret 
the multi-component broad-band spectrum of PG1211-I-143. 
For the highly ionised absorber, the high velocity is again 
constrained by the presence of narrow absorption lines in 
the EPIC data, with the highly ionised XSTAR compo- 
nent having an apparent blue-shift (in the observer frame) 
of ~0.06±0.01, indicating an outflow velocity (in the AGN 
rest frame) of v~0.14±0.01c. For the warm absorber the 
XSTAR fits yield apparent redshifts of ~0±0.02 (pn) and 
~0.02±0.02 (MOS). Constraints on the velocity of the mod- 
erately ionised absorber are less secure, being derived from 
fitting the mid-band spectral curvature and the rather com- 
plex absorption structure below ~2 keV where the ionised 
gas recovers transparency. The column density and ioni- 
sation parameter also affect the spectral curvature while, 
at lower energies, the re-emission spectrum has overlapping 
spectral structure. It does appear, however, that the mod- 
erately ionised gas is outflowing at a somewhat lower, but 
still high velocity. 

In contrast, the emission line spectra in both pn and 
MOS models prefer a low mean velocity in the rest frame of 
the galaxy, with an apparent redshift of the XSTAR emission 
spectra of ~0.07±0.01, corresponding to a mean outflow ve- 
locity in the rest frame of PG1211+143 of ~3000±3000 km 
s~^. Visual examination of the spectral components in fig- 
ures 5 and 6 suggests the velocity sensitivity of the emission 
line gas derives from the strong OVII and OVIII lines at 
-0.5 and ~0.6 keV. 

5.2 Component luminosities, absorption and 
re-emission 

Overall and component luminosities in the pn and MOS 
model fits are sufficiently similar to take their average val- 
ues. For the 2001 observation we find an overall lumi- 
nosity Lo. 4-10-^1. 2 X lO** erg s~^, with a 2-10 keV lu- 
minosity L2-io~5.0 X 10"*^ erg s"'^. The primary power 
law is the dominant component with Lo.4-io~5 x 10*^ erg 



s~^; this increases to —1.3 x 10*'' ergs^^ when corrected 
for absorption (90 % being in the moderately ionised ab- 
sorber) . The secondary power law has an observed luminos- 
ity Lo.4-10— 6 x lO''* erg s~^, increasing to —7 x lO''* erg s~^ 
when corrected for absorption. 

Re-emission from the ionised outflow provides a strong 
contribution to the observed soft X-ray flux below —1 keV, 
with an integrated luminosity in the emission line spectrum 
Lem— 1.8x lO''^ erg s~^ in the moderately ionised component 
and Lem — l.S x 10*^ erg s~^ in the high ionisation compo- 
nent. Comparing the absorbed and re-emitted luminosities 
we estimate the covering factors, for an optically thin ra- 
dial flow, of CF— 0.2 and CF— 0.1 for the warm and highly 
ionised gas components, respectively. 



6 CHECKING FOR CONSISTENCY WITH 
THE RGS SPECTRUM OF 2001 

A key constraint of the new spectral model is to produce 
a strong soft excess with no strong narrow spectral lines in 
absorption or emission, although there should be weak ab- 
sorption lines imposed in the soft band by the highly ionised, 
high velocity flow. Several such lines, at low signal-to-noise 
in the RGS data, were claimed in our initial analysis re- 
ported in P03. 

Examination of figures 5 and 6 suggest the narrow ab- 
sorption line spectrum would be significantly 'diluted' by 
the dominance, at low energies, of the weakly absorbed 
secondary power law continuum. However, the absence of 
strong, narrow emission lines in the RGS spectrum requires 
substantial line broadening of the re-emission component in 
the EPIC model spectrum. Conceivably this could be a re- 
sult of a high velocity flow integrated over a wide angle, 
which would be consistent with the estimated covering fac- 
tor of —0.2, for a partially fllled outflow. 

The high intrinsic resolution of the RGS therefore offers 
a valuable check on the model spectrum of PG1211+143 over 
the soft X-ray band, even if the emission lines are strongly 
broadened. Testing the model against the RGS data from 
the 2001 observation initially gave an unacceptable (x^ per 
dof = 845/636), with the steep power law (F— 3) dominant. 
Examination of the datamodel residuals showed this to be 
due to substantial low energy structure, too broad to be fit- 
ted by the narrow line input spectrum of the XSTAR model. 
We therefore added a gaussian smoothing factor to the emis- 
sion line component in the model (gsmooth in Xspec). The 
broadened emission lines now increased in strength to match 
the observed spectral structure, with the power law compo- 
nent becoming correspondingly weaker, yielding an excellent 
fit (x^ per dof = 649/635 over the 0.35-1.5 keV band). Both 



© 2006 RAS, MNRAS 000, 1-?? 



6 K. A. Pounds et al. 




i! d 



ij-f^. 



0.5 




I 



^ r 



observed energy (keV) 

Figure 7. The spectral model fitted to the 2001 RGS data of 
PG1211+143. The upper panel shows the fit with no added line 
broadening; the lower panel includes a gaussian smoothing of the 
emission line spectrum providing a much better fit as described in 
Section 6. The RGS data and overall model fits are shown in red 
and black, while the separate model components are the absorbed 
primary power law (green), secondary power law (dark blue) and 
input emission line spectra for the warm (light blue) and highly 
ionised (magenta) gas 



unsmoothed and smoothed spectral fits are reproduced in 
figure 7, where the soft X-ray emission can be identified in 
the figure, in order of increasing energy, with K-shell emis- 
sion of N, O, and Ne, together with several strong Fe L lines. 
The apparent redshift from the XSTAR emission spectrum 
fit of ~0.08±0.01 was consistent with the EPIC data fits, 
while - encouragingly - the RGS data also requires a steep 
power law component (rrN^3.3±0.2). 

Figure 8 provides another view of the spectral structure 
in the soft X-ray band, showing the broad spectral lines that 
are successfully modelled by re-emission from the ionised gas 
in PG1211-I-143. The best-fit gaussian (cr~25 eV at 0.6 keV), 
would correspond to velocity broadening of 29000 km s^^ 
FWHM, which again would suggest the high velocity fiow 
occurs over a wide angle. Much better data will be required 
to determine the true emission line profiles and constrain 
the effects of saturation in the principal resonance lines. 

The important point for the present study is that the 
spectral model developed to fit the EPIC data is fully con- 
sistent with the RGS data provided the emission lines are 
strong but broad. In addition, the RGS spectral fit requires 
a steep power law continuum, making a significant contri- 
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Figure 8. Upper panel) Ratio of the RGS data to the model 
described in Section 6. The mid panel shows the data:model ratio 
when the metals C-Mg are removed from the XSTAR emission 
line spectra, revealing 4 broad peaks corresponding (from the left) 
to resonance line emission from NVII, OVII, OVIII and NelX. 
The lower panel shows the same dataimodel ratio when Fe is also 
removed from the emission line spectra, yielding a further broad 
excess at ~0.7 keV due to Fe-L emission 



bution to the observed soft X-ray fiux below ~1 keV. The 
integrated luminosity in the RGS emission line spectrum, of 
Lem^l.S X 10*'' erg s~^, is consistent with the values de- 
rived from the EPIC spectral fits. To place the soft X-ray 
emission in a wider context we note it represents some 0.5% 
of the bolometric luminosity of PG1211-I-143. In compari- 
son, an XMM-Newton observation of NGC1068, a Seyfert 2 
galaxy with a similar bolometric luminosity, found an inte- 
grated soft X-ray emission luminosity of only ~5 x 10*^ erg 
s~^ (Pounds and Vaughan 2006). The much stronger soft 
X-ray line emission from PG121H-143 supports the view of 
a strong centrally-condensed emission region that would be 
hidden from view in a Seyfert 2. 



7 APPLYING THE NEW SPECTRAL MODEL 
TO THE 2004 OBSERVATION OF 
PG1211+143 

As the new spectral model was based on the evidence for 
variable ionised absorption, and re-emission, from a direct 
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observed energy (keV) 

Figure 9. The spectral model fitted to the 2004 pn camera 
data of PG1211+143. The model components are primary power 
law (green), additional power law (dark blue), re-emission from 
the moderately ionised gas (light blue) and re-emission from the 
highly ionised gas (magenta) 



comparison of the 2001 and 2004 data sets it is clearly im- 
portant to show that the 2001 model parameters can be 
varied in a physically reasonable way to also describe the 
2004 spectra. 

To that end we tested the new spectral model with the 
2004 data of PG121 1-1-143, keeping the element abundances 
in XSTAR and ionisation parameters unchanged. Justifica- 
tion for the constant ionisation parameters is based on the 
near-identical hard X-ray fluxes (which we take as a good 
proxy for the ionising luminosity) . Figure 9 illustrates the fit 
and the separate model components for the 2004 pn data, 
which yielded an acceptable P^r dof = 691/625. The 
fit to the 2004 MOS data was statistically excellent, with 
per dof = 350/351, as was that to the 2004 RGS data. 
Comparing the 2001 and 2004 model fits to both EPIC and 
RGS data shows that the spectral change is primarily due 
to a stronger secondary power law component, though the 
'warm' absorption of the primary power law and the soft 
X-ray line emission also appear to be weaker. The column 
density of the highly ionised absorber is a factor ~2 lower, 
consistent with the significantly weaker absorption line at 
~7 keV in 2004. Dilution by a stronger steep power law 
continuum and weaker line emission both contribute to the 
reduced soft X-ray structure evident in the 2004 RGS spec- 
trum (figure 10). 

In summary, we find reasonable changes in the 2001 
spectral model also fit the 2004 data. The net effect of a 
stronger secondary power law, with reduced ionised absorp- 
tion (and re-emission) can satisfactorily match the ratio and 
difference spectra shown in figures 3 and 4. However, given 
the relative complexity of the spectral de-convolution we 
stress that more than 2 spectral 'snapshots' will be required 
to fully resolve and understand the spectral variability. For 
example, while our modelling suggests the reduced soft X- 
ray structure in the 2004 spectrum is a combination of a 
stronger continuum and weaker line emission, we would ex- 
pect further observations to show this to be coincidental. 
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Figure 10. Ratio of the 2004 RGS data to the model shown 
in figure 9 when the metals C-Fe are removed from the XSTAR 
emission line spectra. Although spectral structure is seen, it is at 
a significantly lower level than in the corresponding plot for the 
2001 RGS data (figure 8, lower panel) 



8 DISCUSSION 

We have explored a new model to describe the complex 
broad-band X-ray spectrum of the narrow emission line QSO 
PG1211-I-143. The model structure was guided by a compar- 
ison of the spectral data from two XMM-Newton observa- 
tions of PG121H-143 in 2001 and 2004, which suggested 
a reduction in ionised absorption (and re-emission) in the 
latter case. In addition, the model allows for a second con- 
tinuum component, described by a steeper power law, to 
take account of evidence from other AGN studies that a 
steep power law can be a principal contributor to spectral 
variability. Although more complex than the conventional 
description of the X-ray spectrum of PG1211-I-143, where 
the hard power law and strong 'soft excess' are both mod- 
elled by Comptonisation of soft accretion disc photons, the 
new model is encouragingly self-consistent. Thus, the main 
continuum absorption is modelled by moderate ionisation 
(warm) matter of column density Nn'^-i x 10^^ cm~^ and 
log^~1.4, with re-emission from gas of the same ionisation 
parameter dominating the observed spectral structure be- 
low ~1 keV. Modelling the higher resolution RGS data con- 
firms the contributions of the broadened emission lines and 
the steep power law continuum to the soft X-ray spectrum 
of PG1211-I-143. From a comparison of the absorbed and 
re-emitted luminosities we find covering factors of ~0.2 and 
~0.1, respectively, for the warm and highly ionised outfiows. 
In turn, the mechanical energy in the highly ionised outfiow 
is confirmed to be an order of magnitude larger than the 
luminosity in the secondary power law, supporting our con- 
jecture that shocks in the outfiow could drive that additional 
continuum component. 

Application of the 2001 spectral model to the 2004 data 
suggests the primary change is due to an increase in the steep 
power law continuum rather than a decrease in absorption 
of the primary power law. In fact, given the quality of the 
shorter 2004 observation, it is possible the ionised absorption 
did not change, a possibility consistent with the very simi- 
lar ionising fiuxes indicated by the comparable hard X-ray 
spectra which would have required a change to be mainly 
in column density (or covering factor). A more extended ob- 
servation of PG1211-I-143, exploring spectral variability over 
hours to days, would be very interesting. 
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8.1 The ionised outflow 

The spectral curvature near ~1 keV, previously seen as the 
onset of a strong 'soft excess' in PG121H-143, is attributed 
in our model primarily to absorption and re-emission in a 
substantial column density (Nn^/ew x 10^^ cm~^) of mod- 
erately ionised gas. Although not as well constrained as the 
high ionisation gas, which produces an array of narrow ab- 
sorption lines (Pounds and Page 2006), the warm absorber 
also appears to be outflowing at a high velocity. With the 
present 'snapshot' spectra it is not clear how these absorp- 
tion components might be related; perhaps they simply rep- 
resent different matter densities in a common flow. 

Identifying the strong absorption line observed at ~7 
koV with FoXXV Hea (here and in Pounds and Pago 2006) 
increases the velocity of the highly ionised outflow from 
~0.09±0.01c (P03) to v~0.14±0.01c, with a corresponding 
reduction in the ionisation parameter (and column density) 
from the XSTAR modelling. These changes can then be ap- 
plied to the radiativoly-driven wind model discussed in P03 
and King and Pounds (2003). 

With an ionising X-ray luminosity (>7 keV) of 
3x10''^ erg s~^ and ionisation parameter ^(=L/nr'^)~1000, 
we have nr^ ~3 x lO"" cm . Assuming a spherically sym- 
metric flow, at an outflow velocity of 0.14c, the mass loss 
rate is then of order M~ 356Mo yr~^, where 6 <1 allows for 
the coUimation of the flow. Prom a comparison of absorbed 
and re-emitted luminosities in the present analysis we found 
b~0.1, yielding an outflow mass rate M ~ 3.5M0 yr~^. 
This compares with MEdd ~ I.6M0 yr^^ for a non-rotating 
SMBH of mass ~4 x lO^M© (Kaspi et al. 2000) accreting at 
an efficiency of 10%. 

The mechanical energy is now ~10*^ erg s^^, a factor 
-^15 larger than the luminosity in the steep power law com- 
ponent. The ratio of mechanical energy to the bolometric 
luminosity remains (eis in P03) roughly consistent with the 
value v/c predicted for a radiatively driven outflow (King 
and Pounds 2003). Again, if the higher velocity equates 
to the escape velocity at the launch radius Kiaunch (from 
an optically thick photosphere or radiatively extended in- 
ner disc), v~0.14c corresponds to 'Riaunch ~ 507?s (where 
-Rs ~ 2GM/(? is the Schwarzscliild radius), or 3 x lO^** cm. 
The EPIC data show significant flux variability in the harder 
(2-10 keV) band on timescales of 2-3 hours (fig 1 in P03), 
which would be compatible with the above scale size relating 
to the primary (disc/corona) X-ray emission region. 

Below ~1 keV we find a similar timescale of variability, 
but with reduced amplitude. The simplest interpretation is 
that the photoionised emission is constant over hours, im- 
plying that the secondary power law component is also pro- 
duced close to the inner disc. Bettor data will bo needed 
to clarify the time variability and scale size of the separate 
model components. 

We have very little information on the geometry of the 
warm outflow, responsible for most of the continuum absorp- 
tion (and rc-cmission) , except that it presumably originates 
close to the black hole (from its apparent high velocity) and 
is largely confined within a radius of ~1 pc (by comparison 
with the much weaker soft X-ray omission scon from Soyfort 
2s). The observed soft X-ray luminosity of Lem'~2 x 10*^ erg 

corresponds to an emission measure En^.l/~10®® cm~* 
for a solar abundance gas with ionisation parameter ^~25. 



With an ionising luminosity ~10'''' erg s~^, the ionisation 
parameter gives n.r'^~4 x lO**^. Assuming a spherical shell of 
radius r, thickness 5v and density n, a second constraint is 
the meeisured column density n.Jr ~ 4 x 10^^ cm~^. Interest- 
ingly, these parameters match the above emission measure 
independent of the assumed scale size. For example, assum- 
ing r^lO^'^ cm, we find n~4 x 10^° cm~^ and Sr ~10^^ cm. 
For r~10^^ cm, an alternative set of self-consistent parame- 
ters are n~4 x 10* cm~^ and 5r ~10^* cm. 

Although the above estimates arc very crude, wo note 
the relatively high densities are essential for much of the 
gas to remain in a state of moderate ionisation so close to 
the powerful continuum source. A likely geometry may be 
for the warm gas to exist in many, dispersed small clouds 
rather than the thin spherical shell assumed above. 

That would extend the picture developed in P03 and 
King and Pounds (2003) to incorporate higher density 
'clouds' entrapped in the fast outflow and responsible for 
the strong low energy absorption (and re-emission) previ- 
ously interpreted as the 'soft excess'. A possible scenario 
might be where an inhomogeneous flow accretes through the 
inner disc to a radius R where radiation pressure causes the 
matter to be launched at the local escape velocity. As the 
outflow expands outward the moan density will fall, as will 
the filling factor of the cooler, more opaque, matter. If the 
secondary power law is indeed formed by shocks in the fast 
outflow, then a centrally condensed absorber could explain 
the stronger absorption we find on the primary power law 
continuum (assumed to form in the disc/corona). 

8.2 The second power law component 

The incentive for including a second power law component 
derived mainly from the detailed studies of MCG-6-30-15 
(Vaughan and Fabian 2004) and IH 0419-577 Pounds et al. 
2004), where difference spectrum analyses showed the pri- 
mary spectral change in each case could be attributed to a 
variable flux, steep power law. It is also interesting to re- 
call that a 'broken power law', with a steeper component 
at lower energies has often been used to parameterise AGN 
spectra, without any physical explanation being offered. In 
the present study of PG1211-I-143 we find the inclusion of a 
second continuum component provides the simplest variable 
to allow a good fit to both 2001 and 2004 EPIC data. In 
contrast, a reduction in the covering factor (modelled by a 
lower column) of ionised absorber on the primary power law 
is not able - alone - to fit both data sets. 

A further reason to include the second power law com- 
ponent in our modelling was to quantify the conjecture that 
the the mechanical energy in the high velocity outfiow in 
PG1211-I-143 might power a continuum component addi- 
tional to primary disc/coronal emission. A key factor in de- 
termining the mass and energy in the outflow is the degree 
of coUimation. Our modelling of the broad-band spectrum 
of PG1211+143 has quantified both the absorbed and re- 
emitted fluxes for the ionised outflow and thereby allowed 
the covering factors (or fractional solid angle) of to be esti- 
mated. With a covering factor CF~0.1, we have seen that 
the mechanical energy in the fast, highly ionised outflow is 
ample to power the second continuum (steep power law) 
component. Few hours variability below ~1 keV, although 
of lower amplitude than in the harder X-ray band, suggests 
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the second continuum component originates at a small ra- 
dius. If powered by the fast outflow, possibilities might be 
internal shocks in the flow (perhaps analogous to the pro- 
cess suggested in Gamma Ray Bursts), or by running into 
the slower moving clouds providing the bulk of the contin- 
uum opacity. A possible X-ray emission mechanism in either 
case could be Comptonisation of opticaJ-UV disc photons, 
where the steep power law would be a result of the relatively 
lower energy content of the electron 'corona' compared with 
that responsible for the primary power law continuum. 



9 A REVISED ASSESSMENT OF AGN X-RAY 
SPECTRA 

In contemplating how widely applicable our new model 
might be to describing the X-ray spectra of type 1 AGN, 
we note that a bolometric luminosity of order 4x10*^ erg 

and reverberation mass estimate for the SMBH in 
PG121H-143 of - 4 X 10^ M© (Kaspi et al. 2000) indi- 
cates PG1211-I-143 is accreting close to the Eddington rate. 
Additional support for that view comes from the optical 
classification of PG1211-I-143 as a Narrow Line Seyfert, now 
thought to be a characteristic of a high accretion ratio (eg 
Pounds and Vaughan 2000 and references therein) . King and 
Pounds (2003) provided a simple physical model whereby 
massive, high velocity outflows can be expected in AGN ac- 
creting at or above the Eddington limit. PG1211-I-143 may 
therefore be a special case among bright, nearby AGN. 

However, it is intriguing that the X-ray spectra of 
many bright type 1 AGN have a similar profile to that 
of PG1211-t-143 shown in figure 1, suggesting ionised ab- 
sorption is a common cause of the 'soft excess'. An alterna- 
tive model that works well for weaker soft excesses invokes 
strong photoionised reflection (Crummy et al. 2006), where 
the smooth spectral profile is explained by relativistic broad- 
ening in the inner disc. 

E\irther observations of PG1211-I-143, extended over 
several days, arc essential to test the new spectral model, 
in particulax to resolve the contributions of variable absorp- 
tion (and re-emission) from ionised gas to the overall spec- 
tral form. The variability timescale of the secondary power 
law will constrain the scale of the region where we predict 
the fast outflow undergoes internal shocks. Detecting hard 
X-ray emission from PG1211-I-143 above ~20 keV would also 
support the need for a continuum component with photon 
index less steep than that in the single power law/ absorp- 
tion models (Gierlinski and Done 2004, Schurch and Done 
2006). Finally, much deeper RGS spectra are needed to re- 
solve the broad emission line profiles indicated in the model 
fits. 



10 SUMMARY 

(1) Previous analyses of the 2001 XMM-Newton observation 
of the bright quasar PG1211-fl43 have reported evidence 
of a high velocity ionised outflow (P03, Pounds and Page 
2006)). By deconvolving the absorption and emission lumi- 
nosities in the broad band spectrum, we now flnd a covering 
factor of order ~0.1 for the fast outflow, confirming the mass 



rate and mechanical energy to be comparable to the accre- 
tion rate and bolometric luminosity, respectively 

(2) We show that an additional outflow component of 
less highly ionised (and hence more opaque) gas can im- 
pose significant curvature on the emerging spectrum near 
-^1 keV, thereby replicating a 'soft excess' when applied to 
a continuum steeper than the canonical value. 

(3) Narrow absorption lines in the EPIC spectra show 
the outflow to be approximately radial, in contrast with the 
assumption of a rclativistically smeared absorber in the soft 
excess modelling of Schurch and Done (2006) . 

(4) However, like Shurch and Done (2006) , we find that 
rc-cmission of the absorbed continuum also makes a sig- 
nificant contribution to the observed soft X-ray flux. The 
high integrated luminosity of the photionised outflow in 
PG1211-I-143 is more than an order of magnitude larger than 
that of NGC1068, a Seyfert 2 galaxy of comparable bolomet- 
ric luminosity, consistent with an expanding flow originating 
close to the SMBH. 

(5) Inclusion of a secondary continuum component, 
which energetically could be powered by the high velocity 
outflow, is supported by a fit to the soft X-ray spectrum of 
the RGS. The RGS analysis also requires the soft X-ray line 
emission to be strongly broadened, possibly by a combina- 
tion of velocity broadening and line saturation. 

(6) Testing the resulting model with data from a second 
XMM-Newton observation of PG1211-t-143 in 2004 shows 
that the spectral change is dominated by an increase in the 
secondary power law component, which mimics the effect of 
reduced ionised absorption of the primary continuum. 

(7) A further interesting consequence of the double 
power law model is to remove the conceptual difficulty in 
considering partial covering as an alternative to the extreme 
relativistic Fe K emission line apparent in a simple power law 
fit (P03). 

(8) Wc suggest that the above model could be gener- 
ally applicable to type 1 AGN accreting at or above to the 
Eddington limit. 
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